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Abstract 



CASCADE is a full hadron level Monte Carlo event generator for ep, ^p and pp and 
pp processes, which uses the CCFM evolution equation for the initial state cascade in 
a backward evolution approach supplemented with off - shell matrix elements for the 
hard scattering. A detailed program description is given, with emphasis on parameters 
the user wants to change and common block variables which completely specify the gen- 
erated events. 



PROGRAM SUMMARY 



Title of Program: CASCADE 2.2.0 



Computer for which the program is designed and others on which it is operable: any with stan- 
dard Fortran 77 (g77 or gfortran), tested on SGI, HP-UX, SUN, PC, MAC 



Programming Language used: FORTRAN 77 
High-speed storage required: No 
Separate documentation available: No 

Keywords: QCD, small x, k t - factorisation, CCFM, parton showers, leptoproduction, pho- 
toproduction, pp- and pp-scattering, heavy quark production, unintegrated PDFs. 

Nature of physical problem: High-energy collisions of particles at moderate values of the frac- 
tional momentum x are well described by resummation of leading logarithms of transverse 
momenta (a B In Q 2 ) n , generally referred to as DGLAP physics. At small x leading-logs of lon- 
gitudinal momenta, (a s In x) n , are expected to become equally if not more important (BFKL). 
An appropriate description valid for both small and moderate x is given by the CCFM evolu- 
tion equation, resulting in an unintegrated gluon density A(x, k t , q), which is also a function 
of the evolution scale q. 



Method of solution: Since measurements involve complex cuts and multi-particle final states, 
the ideal tool for any theoretical description of the data is a Monte Carlo event generator 
which embodies small-x resummation, in analogy to event generators which use DGLAP 
resummation. The CCFM evolution equation forms a bridge between the DGLAP and BFKL 
resummation and can be applied to generate the initial state branching processes. The CCFM 
equation can be formulated in a way suitable for carrying out a backward evolution, which 
is an essential requirement to efficiently generate unweighted Monte Carlo events. 

Restrictions on the complexity of the problem: The following hard subprocesses can be simu- 
lated: 7*5* -»• qq{QQ), 75* -> J 1^9, 9*9* -> qqjQQ), 9*9* -> J 1^9, 9*9* -> X, 9*9 -> 99, 
9*q -»• 59/ 5*5* -> h°, g*g -> 5 *5* -> ^Q<2, 5*5* -> ^99, 5*5* -> Wttfc. 

The present version is applicable for HERA, TEVATRON and LHC processes. 

Other Program used: Pythia (version > 6.4) for hadronisation, Bases/Spring 5.1 for inte- 
gration (supplied with the program package). 



Unusual features of the program: None 



1 The CCFM evolution equation 

The formulation of the CCFM HHU parton evolution for the implementation into a full 
hadron level Monte Carlo program is described in detail in (5K6). Here only the main re- 
sults are summarized and discussed. The pattern of QCD initial state radiation in a small- x 
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Figure 1: Kinematic variables for multi-gluon emission in leptoproduction (a) and hadroproduction 
(b). The t-channel gluon four-vectors are given by k n and thegluons emitted in the initial state cascade 
have four-vectors p n . The maximum angle (a function of the rapidity) for any emission is obtained 
from the quark box, as indicated with 3. 

event in ep and pp(pp) collisions is illustrated in Fig. [T] together with labels for the kinemat- 
ics. According to the CCFM evolution equation, the emission of partons during the initial 
cascade is only allowed in an angular ordered region of phase space. In terms of Sudakov 
variables T and 3 the quark pair momentum is written as: 

Pq +Pg = T(p^ +Ep^) + Q t , (1) 

where p^ and p^ are the four-vectors of incoming particles (electron-proton, proton-antiproton 
or proton-proton), respectively and Qt is the transverse momentum of the quark pair in the 
center of mass frame of p^ and p( 2 > (cms). The variable 3 is related to the rapidity Y in the 
center of mass (CMS) frame via 

r _Wf±£0=W>y (2 ) 



2 °\E-p z J 2 

Using E = E q + Eq and p z = p q Z +Pq z gives E + p z = Tyfs, E -p z = TE^ with E = \fsj1 
and s = (p^ + p( 2 )) 2 being the squared center of mass energy. Therefore 3 can be used to 



define the maximum allowed angle in the evolution. The momenta pi of the gluons emitted 
during the initial state cascade are given by (here treated massless): 



Pi = v i (pW+Z i pW)+p ti , & = % (3) 

svf 

with Vi = (1 — Zi)xi-\ and Xi = z%Xi^\. The variables X{ and Vi are the momentum fractions 
of the exchanged and emitted gluons, while zi is the momentum fraction in the branching 
(i — 1) — > i and pu is the transverse momentum of the emitted gluon i. Again the rapidities 
Ui are given by yi = —0.5 log & in the CMS frame. 

The angular ordered region is then specified by (Fig. Hfr and the lower part of the cascade 
in Fig. [It, for the upper part the variables have to be changed accordingly): 

£o < Ci < • • • < £n < S, (4) 

which becomes: 

Zi-iqi-i < Qi, (5) 
where the rescaled transverse momentum qi of the emitted gluon is defined by: 

qi = Xi-i ^fs£i = Pu . (6) 

1 — Zi 

The scale q (related to the maximum angle) can be written as: 

q 2 = ? 2 Zs = s + Ql (7) 

with s = (p q + pq) 2 and the relation of q to a particular choice of the factorisation scale ///in 
the collinear approach becomes obvious. 

The CCFM evolution equation can be written in a differential form Q, which is best 
suited for a backward evolution approach adopted in the Monte Carlo generator Cascade 

SB 

_ 2 d xA(x,kt,q) f j d(f> P(z, q/ z, k t ) 
Q dq 2 A s (q,Q ) 

where A(x,kt,q) is the unintegrated gluon density depending on x, k t and the evolution 

variable q. The splitting variable is z = x/x' and k t = (1 — z)/zq + k t , where the vector q is 
at an azimuthal angle <\>. The Sudakov form factor A s is given by: 

with a s = AOla = For inclusive quantities at leading-logarithmic order the Sudakov 
form factor cancels against the 1/(1 — z) collinear singularity of the splitting function. 

The original splitting function P g (zi,qi, k t i) for branching i is given by (neglecting finite 
terms as they are not obtained in CCFM at the leading infrared accuracy (cf p. 72 in |3J): 



P 3 (zi,mMi) = a,{q \ {1 _ 2,)2) + 2iMzA„fe,,, 2 ,4), (10) 

1 Zi Zi 



where the non-Sudakov form factor A„, is defined as: 



IIS 

f 1 dz' f dn 2 

log A ns = -a s (k 2 ti ) / — / -\Q(k ti - q)G(q - z'q tl ). (11) 

Jo z J q 



The implementation of the full splitting function including non singular terms can lead 
to inconsistencies. Replacing naively only — > — 2 + z(l — z) in the CCFM splitting 
function can lead to negative branching probabilities. 

In (7| it was suggested to use: 

P{z,q,k) = a s (k 2 )(^^ + (l-B)z(l-z)^A ns (z,q,k) (12) 

+a s ((1 - zfq 2 ) + Bz(l - z)\ , 

where B is a parameter to be chosen arbitrarily between and 1, we take B = 0.5. As 
a consequence of the replacement, the Sudakov form factor will change, but also the non- 
Sudakov form factor needs to be replaced by: 

log A ns = -a 8 (kf) £ dz' (^—^- + (1 - B)z(l - z)j J ^-&(k - q')@(q' - zfq). (13) 

2 Backward evolution: CCFM and CASCADE 

The idea of a backward evolution [8,91 is to first generate the hard scattering process with the 
initial parton momenta distributed according to the parton distribution functions. This in- 
volves in general only a fixed number of degrees of freedom, and the hard scattering process 
can be generated quite efficiently. The initial state cascade is generated by going backwards 
from the hard scattering process towards the beam particles. 

According to the CCFM equation the probability of finding a gluon in the proton de- 
pends on three variables, the momentum fraction x, the transverse momentum squared k 2 
of the exchanged gluons and the scale q = x n y/sE, which is related to the maximum angle 3 
allowed for any emission. To solve eq.© the unintegrated parton distribution A(x, k t ,q) has 
to be determined beforehand. 

Given A(x, k t , q), the generation of a full hadronic event is separated into three steps, as 
implemented in the hadron level Monte Carlo program CASCADE: 

• The hard scattering process is generated, 

a = J dkt 1 dkt 2 dxidx2A(xi,kti,q)A(x2, k t2 ,q)d-(ki + k 2 ->■ X) , (14) 

with ki(k 2 ) being the momenta of the incoming partons to the subprocess k± + k% — >• X 
with X being the final state. The definition of a follows [10J . The available processes 
are shown in Tab [21 The momenta of the incoming partons are given in Sudakov rep- 
resentation: 

kt = xip^ + X2V^ + h c± xip^ + kti 
k 2 = Xip^ + X2P^ + h ~ X2P y2 ' ) + k t2 

where the last expression comes from the high energy approximation {x\ t 2^1)/ which 
then gives — k 2 ~ k 2 . 

• The initial state cascade is generated according to CCFM in a backward evolution ap- 
proach. 

• The hadronisation is performed using the Lund string fragmentation implemented in 
Pythia /Jetset fll). 



The parton virtuality enters the hard scattering process and also influences the kinemat- 
ics of the produced particles (Zo, W, Higgs and quarks) and therefore the maximum angle 
allowed for any further emission in the initial state cascade. This virtuality is only known 
after the whole cascade has been generated, since it depends on the history of the parton evo- 
lution (as x in eq.([T5j> may not be neglected for exact kinematics). In the evolution equations 
itself it does not enter, since there only the longitudinal energy fractions z and the trans- 
verse momenta are involved. This problem can only approximately be overcome by using 
k 2 = kf/(l — x) for the virtuality which is correct in the case of no further parton emission in 
the initial state. 

The Monte Carlo program Cascade can be used to generate unweighted full hadron 
level events, including initial state parton evolution according to the CCFM equation and 
the off-shell matrix elements for the hard scattering process. It is applicable for pp, pp, pho- 
toproduction as well as for deep inelastic scattering. A discussion of the phenomenological 
applications of CASCADE can be found in fl2|| . 

The typical time needed to generate one event is similar to the time needed by standard 
Monte Carlo event generators like Pythia )HT| . 

2.1 The unintegrated parton density 
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Table 1: Recommended CCFM unintegrated parton distribution functions included in CASCADE. 
See also sec. 14.5.51 

The CCFM unintegrated parton density xA(x, h, q) can be obtained from a forward evo- 
lution procedure as implemented in SMALLX [1611171 by a fit to the measured structure func- 
tion F2 as described i.e. in USUI. From the initial parton distribution, which includes a Gaus- 
sian intrinsic k t distribution, a set of values x and k t is obtained by evolving up to a given 
scale log q using a forward evolution procedure. Technically the parton density is stored on 
a grid in log x, log k t and log q and a linear interpolation is used to obtain the parton density 
for values in between the grid points. The data file (i.e. ccf m-xxxx . dat) containing the 
grid points is read in at the beginning of the program. 

Several sets (J2003 set 1 - 3 [15] ( IGLU=10 01-10 3), set A Ha IGLU=1010-1013 and 



set B |jl3| IGLU=1020-1023) of unintegrated gluon densities are available with the input 
parameters fitted to describe the structure function F^x, Q 2 ) in the range x < 5 • 1CT 3 and 
Q 2 > 4.5 GeV 2 as measured at HI IMEIl and ZEUS lIMi Set JS E) (IGLU=1) is fitted only 
to F2 (x, Q 2 ) of Ref. [18J. The collinear cutoff kf 11 = Qq which regulates the region of z — > 1, is 
applied both to the real emissions as well as inside the Sudakov form factor. In JS, J2003 set 1 
- 3 and set A we have kf 11 = Qq = 1.3 GeV. Similarly, fits can be obtained using different val- 
ues for the soft cut kf ut = 0.25 GeV, which are available in set B. The set C (I GLU= 1101) flU 
uses the full splitting function, as described in eq. ((T2|) , with a value for ^q CD = 0.13 GeV. 
This set was obtained by fitting simultaneously the inclusive F2 (x, Q 2 ) and jet measurements 
in DIS resulting in a changed intrinsic k t distribution. A CCFM parametrisation of the va- 
lence quark distribution is available (using as starting distribution CTEQ 5 |22| and evolved 
with a splitting function P qq |3) including angular ordering of the emitted gluon). The avail- 
able CCFM uPDF sets with the parameters of the starting distributions are listed in Tab. |2.1[ 

With the parameter IGLU also other unintegrated gluon densities are accessible: a sim- 
ple numerical derivative of a standard integrated gluon density dxg fcP ^ taken from [23J 
(IGLU=2), the one in the approach of Blumlein (241 and coded by |2§J|26l (IGLU=3), the 
unintegrated gluon density of KMS0 p7| (IGLU=4, stored in kms . dat), the one of the satu- 
ration model by |28~| (I GLU= 5, i.e. parameter set including charm jl and the one ofKMFEIl29j 
(IGLU=6, stored in kmr . dat). 

Initial state parton showers can be only generated for the CCFM unintegrated gluon den- 
sity (with IGLU=1 and IGLU=1001-1023, 1101). For all other sets only the cross section 
can be calculated without explicit inclusion of initial state parton showers, since the angu- 
lar variable, essential for angular ordering in the initial state cascade, is not available in the 
uPDFs. However, the transverse momenta of the incoming partons are properly treated. 
Only the KMR set (IGLU=6) provides a prescription for the emission of at most one addi- 
tional gluon. 

3 Hard processes in CASCADE 

Different sets of hard processes applicable for lepto (photo) - and hadroproduction have been 
calculated and are implemented in CASCADE. The available processes are listed in Tab.|2l 

3.1 Lepto(photo)production 

CASCADE can be used to simulate leptoproduction events over the whole Q 2 range. By 
fixing the light quark masses to m q = 0.25 GeV and a s for small fi, the hard scattering matrix 
element remains finite over the full phase space. The total cross section is simulated by 
selecting IPRO=10 and NFLAV=4 (5) . With IPRO=10 light quarks (u,d,s) are selected and 
with NFLAV> 3 the program automatically includes heavy flavour production via the process 
IPRO=ll and IHFLA=4 up to IHFLA=NFLAV. The flag IRE1 indicates, whether beam 1 has 
a internal structure: IRE 1=1 is used to generate resolved photon events. 

Heavy flavour production can be generated separately via IPRO=l 1 . The value of IHFLA 
determines the heavy flavour to be generated. 

The matrix element for jg* — > J/ip(T)g calculated in [31-33] is available for quasi-real 7's 
via the process IPRO=2. The flavour of the Onium is selected via IHFLA, i.e. IHFLA=4 for 

1 A. Stasto kindly provided the program code. 

2 The values of a s and quark masses of [28 1 are not automatically used in the cross section calculation, but 
need to be set explicitly. 

3 M. Kimber kindly provided the program code . 
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Table 2: Processes included in Cascade. Q stands for heavy quarks, qfor light quarks. 

J/ip and IHFLA=5 for T. The matrixelement including J/ip(T) polarisation and subsequent 
leptonic decay can be selected with IPSIPOL=l. 

CASCADE can be used to simulate real photoproduction events by using KBE1=2 2. The 
same options as for leptoproduction are available. Resolved photon events can be generated 
with IRE1 = 1. 

3.2 Hadroproduction 

The hadroproduction processes available are listed in Tab [2] The flavour code for beam 1 (2) 
can be chosen as KBE 1=2 2 1 2 for proton or KBE l=-2 2 1 2 for anti - proton, for beam 2 KBE2 
is changed accordingly. 

CASCADE can be used to simulate heavy quark production in pp or pp collisions (5*5* — > 
QQ IPRO=ll for heavy flavour production, and IHFLA=4 (5) for charm (bottom) quarks), 
but also for light quarks with IPRO=10. The matrix element for 5*5* — > J/tpg calculated 
in Ij3ll - l33| is available via the process IPRO=2. The matrixelement including J/ip polarisa- 
tion and subsequent leptonic decay can be selected with IPSIPOL=l. The process 5*5* — > x 
is available with IPRO=3 including all three x states with appropriate spin and angular mo- 
mentum. The flavour of the Onium is selected via IHFLA, i.e. IHFLA=4 for J/tp(xc) arid 
IHFLA=5 for T(xb)- 

The process 5*5* — > hP with the matrix element calculated in l|34| is available via IPRO=l 02, 
the Higgs mass can be selected via PMAS (25). 

The process 5*5* — > Zqq, calculated in [35,36], is available via IPRO=5 3 for light quarks 
and I PRO =5 4 for the heavy quarks. The flavor of the heavy quark is selected via 
I HFLA=4 (5,6) for charm(bottom, top). The process 5*5* -> Wqiqj, calculated in fl35M36 |, is 
available via IPRO=5 13 for light quarks and IPRO=514 for heavy quarks, where the flavour 
of the heaviest quark is defined by IHFLA=3 , (4,5). Which of the quarks are produced 
depends on the charge of the W which is randomly selected. 



The processes g*g — > gg and g*q — > gq |38| are included via process IPRO=10. The 
individual processes can be selected for g*g* — > qq via IRPA=1, g*g — > gg via IRPB=1 and 
g*q — )• gq via IRPC=1. Note that here one of the partons is treated on-shell. For the quarks 
the unintegrated quark distribution (for valence quarks) is used. 

The process qg* — > Zq, calculated in |35U37B . is available via IPRO=5 01, using the unin- 
tegrated (valence) quark distribution for the on-shell quark. 

3.3 n, and the choice of scales 

The strong coupling a s is calculated via the PYTHIA |[TT| subroutine PYALPS. Maximal and 
minimal number of flavours used in a s are set by MSTU (113) and MSTU (114) , Aqod = 
P ARU (112) with respect to the number of flavours given in MS TU ( 1 1 2 ) and stored in the 
PYTHIA common block C0MM0N/PYDAT1 /. In the initial state cascade according to CCFM, 
the transverse momenta of the i-channel gluons are allowed to perform a random walk for 
small z values and k t can become very small. In the 1/z part of the splitting function we use 
{j, = kt as the scale in a s (fi) and in the 1/(1 — z) part /i = p t is used. In addition we require 
> Qq, resulting in a s (/i > Qq) < 0.6. 

The scale fi, which is used in a s in the hard scattering matrix element, can be changed 
with the parameter IQ2, the default choice is fi 2 = p\. 

The renormalisation scale dependence of the final cross section can be estimated by chang- 
ing the scale used in a s in the off-shell matrix element. Since here we are using the LO a s 
matrix elements, any scale variation will change the cross section. In order to obtain a rea- 
sonable result, the uPDF was fitted to describe F2 by varying the scale fi r . The set A0-,B0- 
correspond to a scale [i r = 0.5pt whereas set A0+,B0+ correspond to a scale {i r = 2p t . 

In order to investigate the uncertainties coming from the specific choice of the evolution 
scale, another definition is applied, relating the factorisation scale only to the quark (or anti- 
quark): Hf = J^- with p t being the transverse momentum of the quark (anti-quark) and 

z = . The set A1,B1 correspond to a scale fif = ■ In the PDF set C J^q CD was fixed to 

Ag^ = 0.13 GeV. 

3.4 Quark masses 

The quark mass for light quarks (u, d, s) is fixed to m q = 0.25 GeV. This, together with the 
treatment of a s at small scales fi, gives also a reasonable total cross-section for photopro- 
duction at HERA energies. The masses for heavy quarks are given by the JETSET / PYTHIA 
defaults (m c = 1.5 GeV, m& = 4.8 GeV) and can be changed according to the PYTHIA pre- 
scription. 

3.5 Initial and final state parton showers 

Initial state parton showers are generated in a backward evolution approach described in 
detail in J5H6J. The initial state parton shower consists only of gluon branchings and is gen- 
erated in an angular ordered region in the laboratory frame. The gluons emitted in the initial 
state can undergo further timelike branchings. The maximum timelike mass m max is calcu- 
lated using the angular constraint. With this mass, the parton which can further undergo 
a timelike branching is boosted to its rest frame with m max but keeping the original energy. 
The timelike branching is performed with the PYTHIA routine P YS HOW. After successful time- 
like branching, the proper mass is associated to the parton and the kinematics are calculated 
appropriately. Gluon radiation from the valence quarks is also included. 

All parameters (like the scale fi in a s , the collinear cut-off kf 11 = Qq) for the initial state 
cascade are fixed from the determination of the unintegrated gluon density. The transverse 



momenta of the partons which enter the hard scattering matrix element are already gener- 
ated in the beginning and are not changed during the whole initial and final state parton 
showering. 

The final state parton shower uses the parton shower routine PYSHOW of PYTHIA with 
the default scale n 2 = 2 • (m 2 ± + m 2 , j_) (IFIN=1), with m 1 ( 2 ) j_ being the transverse mass of 
the hard parton 1(2). Other choices are possible: fi 2 = s (IFIN=2) and fi 2 = 2 • (m 2 + m 2 ) 
(IFIN=3). In addition a scale factor can be applied: SCAFx^ 2 (default: SCAF=1). 

3.6 Remnant treatment 

In CASCADE version 1 the proton remnant was built in subroutine CAREMN, which is a 
slightly modified version of the PYTHIA /LEPTO subroutine PYREMN. No intrinsic transverse 
momentum, in addition to the transverse momentum from the initial state cascade, was in- 
cluded. 

From version 2.0 on the proton remnant can be generated directly via PYTHIA by select- 
ing ILHA=1 (which is now the default). The structure of the event record is then identical 
to the one obtained from a standard PYTHIA run. 

3.7 Hadronization 

In CASCADE version 1 the hadronization was done exclusively by PYTHIA. From version 
2 onwards, events can be written into a file (via switch ILHA=1) according the LHA Ac- 
cord |39| , which can be read by any hadronization program (like PYTHIA or HERWIG ), gen- 
erating the remnants and performing the hadronization. With I LHA= 1 the hadronization is 
performed within PYTHIA. The old CASCADE format is obtained with ILHA=0. Please note, 
that top decays can only be simulated properly within the PYTHIA fragmentation and are 
therefore available only with I LH A= 1 . 

4 Description of the program components 

In CASCADE all variables are declared as Double Precision. The Lund string model is 
used for hadronization as implemented in PYTHIA fTT|. The final state QCD radiation is 
performed via PYSHOW from PYTHIA . The treatment of the proton remnant follows very 
closely the ones in LEPTO [40J for the leptoproduction case and the one in PYTHIA for the 
proton - proton case. However slight modifications were needed to adapt to the cascade 
treatment here. 

The unintegrated gluon density is stored on data files (cc f m . dat , kms . dat , kmr . dat), 
and is read in at the beginning of the program. 

The program has to be compiled and linked together with PYTHIA 6, to ensure that the 
double precision code of Jetset is loaded. 

4.1 Random number generator 

Since the variables are declared as double precision, also a double precision random number 
generator has to be used to avoid any bias. The function DCASRN gives a single random num- 
ber, the function DCASRNV returns an array of length LEN of random numbers. The default 
random number generator is RANLUX (called in DCASRN and DCASRNV) The source code of 
RANLUX (extracted from CERNLIB) is included in the distribution. The user can change this 
to any preferred Double Precision random number generator. 

4.2 Integration and event generation 

The integration of the total cross section and the generation of unweighted events is per- 
formed with the help of BASES /SPRING (4T| , which is included in the distribution package. 



4.3 Program history 

CASCADE 



Version 2.2.00 (Aug 2010) 

time like showering in initial state cascade 

forward DY process added 

J/psi and chi_c (b) production added 

all commons /variables now in double precision (version before cause 
problems for LHC energies and small x) 
obsolete routines from bases removed 
Bug in kinematics for Higgs prod corrected . 



Version 2.1.00 (31. Dec 2009) 

pt cuts are now done all in meof f sh routines on part on s entering ME calc 
W/Z + QQ production added 

QCD jet production with qg->qg and gg - >gg for on she 11 quark ( on she 11 g) included 



CASCADE 

* Version 2 . . 02-beta (Sept 2009) 

* meof f sh includes again correction for k~ 2 = kt "2 (needed for 

* reasonable description of F2 at small x 

* bug corrected for this correction in px, py (generate phi instead of cos phi ) 

* only one library created, including all files (16. Aug 2008) 

* IPR0=15 (qg ->qg ) added using valence quarks (12. Aug 2008) 

* CCFMSTFQ for the unintegrated valence quark distribution added(12.Aug 2008) 

* PTCU is now applied always in the lab frame, no longer in the CM 

* frame, since divergencies could be still there. (12. Aug 2008) 

* auto too Is now also with shared libraries (20 . Mar 2008) 

* updated steering files, also for GENSER.(20. Mar 2008) 

* bug in event corrected: only events with xsec>0 accepted (5.4.08) 

* NOW CERNLIB free version 

* in caps: for ILHA=10 MSTJ(41)=1 enforced to avoid isolated photons 

* from FPS in pythia event record 

* Simple MC integration and generation routine included via IINT=1 

* Valence quark distributions and g"*q -> gq and g"*g -> gg processes 
are included. 



CASCADE 



Version 2.0.01 (24. Dec 2007) 

improvements done on installation using autotools 
working now with make and make install 



Versi 
LHA i. 
ILHA : 

updat 
updat 
updat 



updat 
meof f 



2.0.00 

nterface for PYTHIA/HERWIG included 

10 uses PYTHIA for final state PS and remnant treatment 
es in caupinit, caupevnt 

in p_semih . F : removed line with P(2,4) =abs (p ( 2 , 3 ) 
in caps.F restored event record also for LST(21)=55 
caused energy-mom mismatch before, 
for ILHA> 1 also use caremn . F 
more precise energy -mom check 
in caremn . F for ILHA> 1 set pt diquark=0 
for ILHA>1 limit chi<0.8 
sh: scale PT2Q changed to average of 2 outgoing quarks 
reduces xsection by ca 20% (compared to vers 1.2) 



Version 12010 

- insq set to 0, as for updf s . . . (in casbran.F) 

- QlamO set to QgO, as in updf .. (in casbran.F) 

- changed upper limit of xsi_hard to min (xsi_hard ( 2 ) , 5 . d8 ) in cascps.F 

showed up in wrong Qmax distribution for PS. 

- no IPS for QG21 (QG22) < Q0 (in cascps.F) 

- final state PS was always essentially switched off by IFINAL=0 in 
caps.F. changed to default IFINAL=1 . 

- ordering in casbran . F was wrong, always only q ordering instead of 
angular ordering (using now casbran-v24 . F) . 



Version 12009 

meof f sh updated: now 2 scales for pp 

unnecessary cut in caps removed (xgtest) , which caused asymmetry 

in part on showering 
bug in cascps . F corrected : xsi_hard was wrong for beam 1 



CASCADE with LHA interface 

* Version 12008 

* ILHA = 1 added: event record written to file in LHA format to be 

* read in by PYTHIA or HERWIG 

* ILHA = 10 added: PYTHIA (using LHA interface) is called directly to 

* do final state PS and hadronsiat ion . 



CASCADE 

* Version 12007 

* date: 2004/12/23 

* bug in cascade for new gluon dat files corrected: 

* IPS was switched off 



bug in steer corrected Nmax now at 1000 



CASCADE (with PYTHIA6.226) 

* Version 12003 

* bug in cascps corrected 

* date: 2004/11/09 13:42:21; 

* steering files updated to new frag . parameters 

* date: 2004/09/10 17:20:26; 

* updated to read scale parameters from gluon file 

* date: 2004/11/08 06:47:41; 

* bugs in cascps and casbran for ppbar corrected 

* bug in casbran (neg t in log) corrected ( found by 



Eduardo Rodrigues ) 



Version 12000 

Higgs production included 

e+e- option and resolved photons included 



Version 10000 

published version for ep and pp 



4.4 Subroutines and functions 

The source code of Cascade and this manual can be found under: 

http : //www .desy.de/" jung/ cascade/ 

C AMAIN main program. 

CAS INI to initialises the program. 

to perform integration of the cross section. This routine has to be called 
before event generation can start, 
to print the cross section and the number of events. 
CAUNIGLU (KF, X, KT, P, XPQ) to extract the unintegrated gluon density xA(x, k t ,q) for a 
proton with KF=22 12, as a function of x =X, fef =KT and q =P. The gluon 
density is returned in XPQ ( ) , where XPQ is an array with XPQ (-6:6). 
to perform the event generation, 
to give a s (n) with \i =RQ. 
to give initial particle and parton momenta, 
to call routines for selected processes: XSEC1. 
to cut on p t for 2—^2 process in integration and event generation, 
matrix element for 7*5* — > qq and q*g* — > qq including quark masses, q 
can be light or heavy quarks, 
matrix element for 7*5* h°. 
four-vector dot product of A and B. 

four- vector dot product of vectors I and J in PY JETS common, 
to generate momenta of final partons in a 2 — > 2 subprocess according to 
phase space. 

to generate kinematics and the event record for ep, ^p and pp processes. 
CAREMN ( IPU1 , IPU2 ) to generate the beam remnants. Copied from LEPTO 6.1 [40] and 

updated for the use in CASCADE. 
CASPLI (KF , KPA, KFSP , KFCH ) to give the spectator KFSP and KFCH partons when a par- 
ton KPA is removed from particle KF. Copied from LEPTO 6.1 [40J and 
updated for the use CASCADE. 
CAPS to generate color flow for all processes and prepare for initial and final 

state parton showers. 
CASCPS ( IPU1 , IPU2 ) to generate initial state radiation. 

COLORFLOW to generate color configuration for g*g — > gg and g*q — > qg processes. 
GADAP Gaussian integration routine for 1-dim and 2-dim integration. Copied 

from LEPTO 6.1 BUI 



CASCADE 



CAEND 



EVENT 
ALPHAS (RQ) 
PARTI 
FXN1 

CUTG (IPRO) 
MEOFFSH 

MEHIGGS 
DOT (A, B) 
DOT1 (I, J) 
PHASE 

P.SEMIH 



4.5 Parameter switches 



BASES/SPRING integration procedure. 
NCAL : (D:=20000) Nr of calls per iteration for bases. 

ACC1 : (D:=l) relative precision (in %) for grid optimisation. 

ACC2 : (D:=0.5) relative precision (in %) for integration. 

Event record output. 
ILHA : (D: = 0) output in LHA accord [393 format 

= 0: CASCADE type output of event record 

= 1: output of event record to be read in by fragmentation programs 

= 10: use LHA format to produce remnant and fragmentation in PYTHIA 

style. 

4.5.1 Parameters for kinematics 

PBEl : (D:=-30) momentum p [GeV/c] of incoming hadron 1 (/INPU/). 

KBE1 : Lund flavour code of incoming hadron 1 (KBE1=11 for electrons, KBE 1=2 2 

for photons, KBE 1 =2 2 1 2 for protons ) 
IRE1 : hadron/lepton 1 has a structure (IRE1=1) or interacts directly with the 

target (IRE 1=0 for a DIS electron) 
PBE2 : (D:=820 momentum p [GeV/c] of incoming proton (/INPU/). 

KBE 2 : Lundflavour code of incoming hadron 2 (KBE 2= 11 for electrons, KBE 2 =2 2 

for photons, KBE2=2 2 1 2 for protons ) 
IRE2 : hadron/lepton 2 has a structure (IRE2=1) or interacts directly with the 

target (IRE2=0 for a DIS electron) 
NFLAV (D: = 5) number of active flavours, can be set by user (/CALUCO/). 

4.5.2 Parameters specific for Ieptoproduction 

QMI : (D: = 5.0) (/VALUES /) minimum Q 2 to be generated 

QMA : (D: = 10 8 ) (/VALUES/) maximum Q 2 to be generated. 

YMI : (D: = 0.0) (/VALUES/) minimum y to be generated. 

YMA : (D: = 1.0) (/VALUES /) maximum y to be generated. 

THEMA, THEMI (D: THEMA = 180., THEMI = 0) maximum and minimum scattering angle 9 
of the electron (/CAELEC/). 

4.5.3 Parameters for hard subprocess selection 

I PRO : (D: = 10) (/CAPAR1 /) selects hard subprocess to be generated. 

=2: 7 g* -> J/1>(T)g, g*g* -> J/^(T)g. 
=3: g*g* -> Xc(b)- 

=10: 7*5* — > qq, g*g* — > qq, g*g — > gg and g*q — > gq for light quarks. 

=11: 7*5* —j- QQ or g * g* -> QQ for heavy quarks. 

=102: g*g* — > h° for Higgs production in hadron-hadron collisions. 

=501: qg* — > Zq for Z+jet production in hadron-hadron collisions. 

=503: g*g* — > Zqqior Z+jet production in hadron-hadron collisions. =504: 

g*g* — > ZQQ The flavour index of the heaviest quark is selected via IHFLA 

=513: g * g* — > Wqiqj for M^+jet production in hadron-hadron collisions. 

=514: g * g* — > WqiQj for VF+jet production in hadron-hadron collisions. 
IRPA: = 1 g*g* — > qq switched on for IPRO=10. 

IRPB: = 1 g*g — > gg switched on for IPRO=10. 

IRPC : = 1 g*q — > gq switched on for IPRO=10. 



IHFLA: 



IPSIPOL: 



PT2CUT (IPRO) 



(D:=4) = 4 flavor of heavy quark produced (in IPRO=10, IPRO=50 4 and 
IPRO=514 ). 

(D:=0) = 1 use matrixelement including J ftp (T) polarisation and subse- 
quent leptonic decay for IPR0=2. 
(D=0.0) minimum p]_ for process IPRO (/CAPTCUT/). 



4.5.4 Parameters for parton shower and fragmentation 



NFRAG : 



IFPS : 



ITIMSHR: 



ICCFM: 



IFIN 



SCAF 



(D: = 1) switch for fragmentation (/CAINPU/). 
= 0: off 
= 1: on 

(D: = 3) switch for parton shower (/CAINPU/). 
= 0: off 

= 1: initial state 

= 2: final state 

= 3: initial and final state 

(D: =1) 

=0: no shower of time like partons 
=1: time like partons may shower 
(D: =1) 

=1: CCFM evolution (all loops) 

=0: DGLAP type evolution (one loop) 

(D:=l) scale switch for final state parton shower 



= 1: /j 2 
= 2: /j? 
= 3: p 2 

(D:=l.) scale factor for final state parton shower 



2(m( t + mi t ) 
s 

2{m\ + m 2 ,) 



4.5.5 Parameters for structure functions, a, and scales 



IRUNAEM: 



I RUN A: 



,(Q 2 



IQ2 : 



(D: = 0) (/CAPAR1/) select running of a e 
=0: no running of a em (Q 2 ) 
=1: running of a em (Q 2 ) 
(D: = 1) switch for running a s . 
=0: fixed a s = 0.3 
=1: running a s (p 2 ) 
(D: = 3) select scale p? for a s (p 2 ). 
=1: p? = 4 • m 2 (use only for heavy quarks!) 



=2: f 
=3: p 2 
=4: p 2 
=5: ix 2 



s (use only for heavy quarks!) 

4 • m 2 + p 2 ± 

Q 2 

Q 2 + p\+ 4 • m 2 



IGLU: 



(D: = 1010) select unintegrated gluon density (/ GLUON/). 

Note that initial state parton showers not possible for I GLU= 2 , 3 , 4 , 

=1: CCFM old set JS2001 

= 1001: CCFM J2003 set 1 CHI 

= 1002: CCFM J2003 set 2 CHI 

= 1003: CCFM J2003 set 3 CHI 

= 1010: CCFM set AO CU 

= 1011: CCFM setA0+d 



= 1012: CCFM set AO- d 

= 1013: CCFM set Al Q3 

= 1020: CCFM set BO d 

= 1021: CCFM set B0+ $M 

= 1022: CCFM set BO- d 

= 1023: CCFM set Bl 

= 1101: CCFM setC CI 

=2: derivative of GRV (23 dx %¥ ] . 

=3: approach of Blumlein [24J 

=4: KMSH2ZI (kms.dat) 

=5: saturation model [28] 

=6: KMR 1291 (kmr . dat) 

4.5.6 Accessing information 



AVGI 
SD 

sss 

PBEAM 
KBEAM 
Q2 
YY 

YY_BAR 
XG 

XG_BAR 
KT2_1, KT2_2 

YMAX, YMIN 
Q2MAX, Q2MIN 
XMAX, XMIN 
AM ( 1 8 ) 
SHAT 

NIA1, NIA2 
NF1, NF2 

Q2Q 

ALPHS 

ALPH 

NIN 

NOUT 



integrated cross section (/CAEFFIC/). 

standard deviation of integrated cross section (/CAEFFIC/). 

squared center of mass energy s (/CAPARTON/). 
energy momentum vector of beam particles (/C ABEAM/), 
flavour code of beam particles (/CABEAM/). 
in leptoproduction: actual Q 2 of 7 (/ CAPAR4 /). 

negative light-cone momentum fraction of parton 1 (7*, g*) (/CASGKI /). 
positive light-cone momentum fraction parton 1 (7*, g*) (/CASGKI/). 
positive light-cone momentum fraction of parton 2 (g*) (/CASGKI /). 
negative light-cone momentum fraction of parton 2 (g*) (/CASGKI /). 
transverse momenta squared k 2 ^ [GeV 2 ] of partons 1(2) which enter to 
the matrix element. 

actual upper and lower limits for y =YY (/ CAPAR5 /). 

actual upper and lower limits for Q 2 (corresponding to KT2_1) of 7 (/ CAPAR5/). 
upper and lower limits for x (/ CAPAR5 /). 

vector of masses of final state particles of hard interaction (/ CAPAR3/). 
invariant mass s [GeV 2 ] of hard subprocess (/ CAPAR5 /). 
position of partons in hard interaction in PY JETS event record (/CAHARD/). 
first and last position final partons / particles of hard interaction in P Y JE T S 

(/CAHARD/). 

hard scattering scale /j, 2 used in a s and structure functions (/ CAPAR4 /). 
actual a s (/ CAPAR2 /). 
a em (/ CAPAR2 /). 

number of trials for event generation (/CAEFFIC/). 
number of successful generated events (/CAEFFIC/). 



4.6 List of COMMON blocks 

COMMON/ CABEAM/PBEAM (2, 5) , KBEAM ( 2 , 5) , KINT (2,5) 
COMMON/ CAHARD/NIA1 , NIA2 , NIR2 , NF1 , NF2 
COMMON/ CAHFLAV/ I HFLA 

COMMON/CAINPU/PLEPIN,PPIN,NFRAG, ILEPTO, IFPS, IHF, INTER, ISEMIH 
COMMON/ CALUCO/KE, KP, KEB, KPH, KGL, KPA, NFLAV 
COMMON/CAEFFIC/AVGI, SD, NIN, NOUT 
COMMON/CAELEC/THEMA, THEMI 
COMMON/ CAGLUON/ I GLU 

COMMON/CAP AR1/IPRO, I RUN A, IQ2, IRUNAEM 
COMMON/ CAP AR2/ALPHS, PI, ALPH, IWEI 
COMMON/ CAP AR3/ AM (18) ,PCM(4, 18) 
COMMON/CAP AR4 /Q2, Q2Q 

COMMON/ CAP AR5/SHAT, YMAX, YMIN, Q2MAX, Q2MIN, XMAX, XMIN 
COMMON/CAP AR6/LST (30) , IRES (2) 
COMMON/CAP ARTON/SSS, CM ( 4 ) , DBCMS (4) 
COMMON/ CAPTCUT/PT2CUT (20) 

COMMON/CASKIN/YY, YY_BAR, XG, XG_BAR, KT2_1 , KT2_2 , PT2H, SHH 
COMMON/ VALUES /QMI, YMI, QMA, YMA 



5 Example Program 

PROGRAM CASMAIN 
Implicit None 
Integer N1,N2 

DOUBLE PRECISION PLEPIN,PPIN 

INTEGER KE, KP , KEB, KPH, KGL, KPA, NFRAG, ILEPTO, IFPS, IHF 
INTEGER INTER, ISEMIH 

INTEGER NIA1, NIR1, NIA2, NIR2, NF1, NF2, NFT, NFLAV 
COMMON/CALUCO/KE, KP, KEB, KPH, KGL, KPA, NFLAV 

COMMON/CAINPU/PLEPIN, PP IN, NFRAG, ILEPTO, IFPS, IHF, INTER, ISEMIH 
C0MM0N/CAHARD/NIA1,NIA2,NIR2,NF1,NF2 
INTEGER I HFLA 
COMMON/ CAHFLAV/IHFLA 

DOUBLE PRECISION THEMA, THEMI , PT2CUT 
INTEGER IRUNA, IQ2, IRUNAEM 
INTEGER IPRO 

COMMON/CAP AR1/IPRO, IRUNA, IQ2, IRUNAEM 

COMMON/CAELEC/ THEMA, THEMI 

COMMON/ CAPTCUT/PT2CUT (20) 

REAL ULALP S , ULALEM 

EXTERNAL ULALP S , ULALEM 

DOUBLE PRECISION QMI , YMI , QMA, YMA 

COMMON/VALUES/QMI, YMI, QMA, YMA 

Integer Iglu 
Common/CAGLUON/ Iglu 

Integer ISEED, I 



I SEED = 124567 
nl = 
n2 = 

C initialize random number generator 
CALL RM48IN (ISEED, N1,N2) 

C initialize PYTHIA 6 parameters 
CALL GPYINI 

C initialize CASCADE parameters 
CALL CASINI 



C Select parton shower (IPS=1 initial, =2 final, 3 init ial+f inal PS ) 

IFPS = 3 
C scale for alpha_s 

C IQ2 =1 mu" 2 = 4m_q" 2 (m_q = light quark or heavy quark depending on IPRO) 
C IQ2 =2 mu~2 = shat 

C IQ2 =3 mu" 2 = 4m_q" 2 + pt "1 (m_q = light quark or heavy quark depending on IPRO) 
C IQ2 =4 mu"2 = q"2 (q"2 of virtual photon) 

C IQ2 =5 mu"2 = q"2 + pt"2 + 4m_q"2 (q"2 of virtual photon) 



IQ2 = 3 

C select process ( IPRO=l for light quarks , IPR0=11 for heavy quarks ) 
IPRO= 10 

C total number of flavours involved 
NFLAV = 4 

C select unintegrated gluon density (D=l) 
Iglu = 1 

C minimum Q"2 of electron to be generated 

QMI = 0.5d0 
C maximum Q"2 of electron to be generated 

QMA = 10D8 
C minimum y of electron to be generated 

YMI=0 . OdO 

C minimum y of electron to be generated 
YMA=1 . OdO 

C maximum theta angle of scattered electron 

THEMA = 180 . ODO 
C minimum theta angle of scattered electron 

THEMI = 0.0D0 
C momentum of beam 1 (electron, proton, antiproton) 

PLEPIN =-27 . 5 

C Lund flavour code for beam 1 ( elect ron=l 1 , photon=22 , proton=22 12 , antiproton=-22 12 ) 
KE=11 

C momentum of beam 2 (proton) 

PPIN = 820. 
C perform fragmentation NFRAG=0/1 

NFRAG = 1 

c for I PRO = 11 which flavour is produced 
IHFLA = 4 

c 

c Start integration of x-section 
c 

CALL CASCADE 

C 

c Print out result of integration of x-section 
c 

CALL CAEND (1) 

c 

c Start event loop 
c 

Do 1=1,100 
c generate an event 
CALL EVENT 

Enddo 

c 

c Print out of generated events summary 
c 

CALL CAEND (20) 

STOP 
END 



6 Program Installation 

CASCADE now follows the standard AUTOMAKE convention. To install the program, do 
the following 

1) Get the source 

tar xvfz cascade-XXXX . tar . gz 
cd cascade-XXXX 

2) Set environment variables for PYTHIA 

example (Please change to the proper path of the libraries) : 
in csh: 

setenv PYTHIA " /home/ jung/cvs/pythia6422 " 
in zsh: 

export PYTHIA= " /home/ jung/cvs/pythia6 422 " 
2) Generate the Makefiles 

./configure — pref ix=install=path — disable-shared 



3) Compile the binary 
make 

4) Install the executable and PDF files 
make install 

4) The executable is in bin 

set the path for the updf data files, if different from 
the default (for example) 

export PDFPATH=/ Users/ jung/ jung/cvs/ cascade2 /cascade-2 . 2 . /share 
run it with: 

cascade < steer_pp-bottom 
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